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Abstract 

A Weyl solution describing two Schwarzschild black holes is considered. We focus on the 
Z2 invariant solution, with ADM mass Madm — "^Mk, where Mk is the Komar mass of each 
black hole. For this solution the set of fixed points of the discrete symmetry is a totally geodesic 
sub-manifold. The existence and radii of circular photon orbits in this sub-manifold are studied, 
as functions of the distance 2L between the two black holes. For L ^ there are two such 
orbits, corresponding to r = 3AIadm and r = 2Madm in Schwarzschild coordinates. As the 
distance increases, it is shown that the two photon orbits approach one another and merge 
when Mk — ^L, where Lp is the golden ratio. Beyond this distance there exist no circular 
photon orbits. The two null orbits delimit a forbidden band for time-like circular orbits, which 
is interpreted in terms of optical geometry. For large L, time-like circular orbits are allowed 
everywhere, as in the analogous Newtonian problem. The analysis is generalised by considering 
a Z2 invariant Weyl solution with an array of N black holes and also by charging the black 
holes, which connects the Weyl solution to a Majumdar-Papapetrou spacetime. 



1 Introduction 



Euler's 3-body problem is a soluble special case of the general 3-body problem, in which a test 
particle moves in the gravitational field of two gravitational sources - point masses - which are 
fixed in space. It is Liouville integrable since, besides the energy and one component of the angular 
momentum, there is a third constant of the motion (see e.g. [H [2]; see [3] for a hyperbolic space 
version of the problem which is still integrable) . A natural relativistic generalisation of this problem 
is to consider the motion of test particles in a Weyl solution describing two Schwarzschild black 
holes, which are kept in equilibrium due to a strut in between them. For this problem there is no 
known analogue to Euler's constant (see [1] for a discussion about this point); equivalently, there is 
no known separability for the geodesic equations. Generic geodesies have, therefore, to be studied 
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numerically. For the special case in which the two black holes have equal mass, there is, however, 
a 1+2 dimensional totally geodesic sub-manifold in which geodesies are Liouville integrable. In 
particular it is straightforward to derive all circular causal orbits. 

Regular (on and outside an event horizon) static black holes in vacuum general relativity do 
not admit causal circular geodesies below a certain critical radius, which lies outside the event 
horizon. Below this radius, the angular momentum term in the radial equation of motion, usually 
centrifugal (with respect to the black hole), becomes centripetal, making equilibrium impossible. 
This change in character of the angular momentum term has a geometric interpretation in terms 
of optical geometry [S]: define "outwards" as the direction in which fixed points sets of the optical 
geometry's isometry group (spheres) increase their proper size; then, the angular momentum term 
is always "outwards" . Similar considerations may be made about stationary black holes. However, 
in this case the optical geometry is a Randers-Finsler rather than Riemannian geometry, making 
the interpretation more subtle [6]. 

In this paper we show that the optical geometry of the aforementioned relativistic Euler's 3- 
body problem is richer than the standard one of the Schwarzschild black hole. To be concrete 
let the two black holes in the Weyl solution considered have equal mass Mk and separation 2L 
between them. This choice of masses implies that there is a discrete Z2 symmetry of the solution. 
The set of fixed points of this discrete symmetry is a totally geodesic sub-manifold, N . The optical 
geometry of this sub-manifold then has the following feature: for sufficiently small but non-zero 
L the notion of "outwards" , defined as above, changes continuously twice, rather than just once. 
This means that, on AA, there is an annulus-like region in which "outwards" means a decreasing 
radial Weyl canonical coordinate. The boundaries of this annular region are the loci of two circular 
photon orbits, and the region itself is a forbidden band for time-like circular geodesies, which are 
allowed everywhere outside this band. As L increases, the two boundaries of the annular region 
approach one another and coalesce when Mk = fL, where 99 is the golden ratio! For larger L, M 
admits no circular photon orbits and time-like circular geodesies are allowed everywhere, just as in 
the Z2 symmetric Euler's 3-body problem. 

Considering a Weyl solution with an array of A^, rather than two, black holes along a line, still 
possessing a Z2 discrete symmetry, leads to two cases. For N even (odd), we essentially recover 
the optical geometry of the N = 2 (N = 1) case. We shall also consider charged black holes in 
Einstein-Maxwell theory, by using a solution generating technique [10] that allows the introduction 
of a charge parameter in Weyl solutions. For black holes, the charge parameter is the ratio of 
charge to mass, which is the same for all the black holes in the solution. In this case, as we increase 
the charge for fixed L,Mk, the proper area of the forbidden band increases; for L = 0, this area 
diverges as we reach extremality, as a consequence of the infinite 'throat' developed by extremal 
black holes. Keeping the charge and Mk fixed, the area of the forbidden band decreases as L is 
increased, just as for the uncharged case, becoming zero at some maximum value of L. 

This paper is organised as follows. After a brief discussion of the Newtonian Euler's 3-body 
problem in section 2, we describe in section 3.1 the Weyl solutions that shall be used in the 
relativistic version of the problem, as well as the circular null and time-like orbits. In section 3.2 
the charged case is considered. In section 4 the interpretation of the forbidden band for time- like 
circular orbits is given in terms of optical geometry. In section 5 we briefly discuss the case of 
multiple black holes. We close with some final remarks. 
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2 The Newtonian version 



In Euler's 3-body problem, a.k.a. the two-centre Kepler problem (see e.g. [7]), with both masses 
equal to M, the motion of a test particle moving in the symmetry plane between the two masses 
obeys 

'dpV . . J'^ 4M 

IE 



2E - V{p) , V{p) 



(1) 



p2 TTTX^ ' 

where p is a radial cylindrical coordinate, E, J are the energy and angular momentum per unit 
mass and L is the distance from either mass to the symmetry plane. We shall use geometrised units 
throughout. Circular orbits are seen by extremising the potential. They obey 



2M 



(p2 + L2)3/2 • 



(2) 



This equation has solution for any p, given L and M. The angular momentum J goes to zero 
(infinity) as p — > (p — > oo) - Fig. [TJ Thus circular orbits exist for any radial distance in 
the symmetry plane. This is therefore the behaviour we expect in the relativistic version of the 
problem for large L/M, and indeed it is the behaviour we shall find. The behaviour for small L/M 
will, however, be quite different. 
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Figure 1: Radial distance, p, of circular orbits in the symmetry plane of the two-centre Kepler 
problem, as a function of their angular momentum, J. The two equal masses, which we set to 
unity, have separation 2L, and various values of L have been used. Observe that circular orbits 
exist for any p. 

In the relativistic problem that shall be considered in the next section, the Newtonian gravi- 
tational potential of one or multiple rods (rather than point-like masses) of mass M, infinitesimal 
width and mass density g plays an important role, since it is a key ingredient in the construction of 
Weyl solutions describing one or multiple black holes. Let us remark that, as a Newtonian problem. 
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the passage from point particles to two equal rods, aligned with the direction of separation, does 
not alter the previous conclusion: in the symmetry plane, circular orbits still exist at any radial 
distance from the symmetry axis. 



3 The relativistic version 
3.1 Uncharged case 

Weyl geometries [HI [9] in four space-time dimensions admit two commuting, mutually orthogonal 
Killing vector fields. In vacuum, the metric may always be written in the form 

The vacuum Einstein equations, R^y = 0, reduce to a harmonic equation in an auxiliary Euclidean 
three space, E^, in cylindrical coordinates (/o, z,^), 

AE3f/ = 0, (4) 

and two partial differential equations for k{p, z) 

d,k = 2pdpUd,u , dpk = p[{dpUf - {d,uf] , (5) 

which become line integrals once the solution for U{p, z) is known. 

A solution with two black holes is obtained by taking U {p, z) to be the potential of two finite 
rods of zero width and linear density q = 1/20 The rods are placed at /> = and in the intervals 
z G [01,02] and z G [03,04]. Then 

p2U{p,z) ^ (-^1 - Cl)(-R3 - C3) ,r.^ 
{R2-Q2){Ra-U)' 

where we have introduced the quantities 



Rk = \Jp'^ + Cl , Ck = z-ak . (7) 

Introducing further 

Yi, = RiRj + CiCj + p" , (8) 

the solution of ([5]), given ([6]), is 

^2k{p,z) ^ ^43^2inil32 _ 

4I42 ^3 1 -Ri -R2 -R3 -R4 

An integration constant could be added to k. We have chosen this constant to be zero. With this 
choice, identifying the azimuthal coordinate with standard period (p ~ (/> + 27r, the metric is smooth, 
on the symmetry axis for z < ai and 2: > 04. For p = and 02 < -z < 03 there is a conical excess, 
given by 

5 = 27r (e-HP=^^-^<^<-^) - l) = 2^ («4 - «3)(a2 - oQ _ ^^^^ 

V / (04 - Oi)(03 - 02) 



^This particular choice is required to avoid curvature singularities at the rod positions. 
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The vacuum solution described by ([3]), ([6]) and ([9|) is actually a three parameters family of 
solutions. Physically, the three parameters may be taken to be the two black holes masses and the 
distance between them. The black hole masses may be computed as Komar integrals: 

= i i = ^ . = ii*-^ = °-^' 

where is the dual 1-form to the time-like Killing vector field d/dt and the 2-surfaces Si correspond 
to the location of the two finite rods in Weyl coordinates. An asymptotic expansion shows that the 
ADM mass is 

Madm = Ml + M2 . (12) 
For the distance we shall take the coordinate distance in Weyl coordinates 

2L = 03 - 02 . (13) 

This is a monotonic function of the proper distance 

ra.3 

Az= ^/gTzdz , (14) 
J 02 

and hence a good measure thereof. 

In order to have a totally geodesic sub-manifold we require Mi = M2 = Mk = Madm/'^- We 
also choose a symmetric coordinate system: 02 = —L and 03 = L. The proper distance between 
the two black holes may then be expressed as 

^ SLjMADM + L)^ f L \ 

(Madm + 2LY + Madm) ' ^ ^ 

where E[x) is a complete elliptic integral of the second kind. It follows that z = is a totally 
geodesic sub-manifold, which has induced metric given by ^ with z = 0. Geodesies in this sub- 
manifold obey 

where the dot represents derivative with respect to an affine parameter and we have introduced the 
energy E and angular momentum J for a particle of mass m. 

3.1.1 Circular photon orbits 

From ()16p . circular null geodesies are determined by the extrema of the potential 

_ e^^(^-°) _ 1 / L+yi^ y ^^^^ 

\Madm + L + ^J{Madm + Lf + p^] 

For L = 0, i.e the single black hole limit, making the coordinate transformation 



, 2Madm 

P = r\n- , (18) 
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the potential becomes 

1 / 2Madm \ 



^(0 = ^[^ ^ ) • (19) 



r 



This is the effective potential for photons in a Schwarzschild black hole and in Schwarzschild 
coordinates. It has the well known extremum at r = 3Madm- Another solution for a constant r 
photon orbit is r = 2Madm, corresponding to the null geodesic generator of the horizon. Thus, in 
Weyl coordinates we have two extrema, for L = 0, at 

PCNO^ (L = 0) = , PCNO2 {L = 0) = VsMadm ■ (20) 

For L ^> Madm, yip) has no extrema, since it is approximately V{p) ~ The interpolation 

between these two behaviours is as follows - Fig[2j As L grows, the two extrema (j20p approach one 
another, in both coordinate and proper distance, and merge for L = Lmax- Beyond this value of 
L, V{p) has no extrema. To see the value of L^ax^ observe that the extrema of (jl7p are given by 

f{p, Madm, L) = - , f{p, Madm, L) = ^^-^dm + L _ L ^ 

To solve for p it is convenient to introduce = p^ / {Madm + -^)^ and y = Madm/L. Then the 
last equation becomes 

111 

- = , — . (22) 

2 ^/lT^ v^l + x2(l+y)2 ^ ' 

Take (j22|) to define a curve y = y{x), between x = and x = \/3; at these values of x, y — > 00. 
These two "points" of the curve correspond to pcNOi and PCNO2 given in (f20l) . Constant y slices of 
the curve have two solutions for y > ymin, one solution for y — ymin and no solution for y < ymin- 
To determine ymin, it is convenient to introduce 

Vminf^^ . (23) 

Then, differentiating (|22p . extremising y and replacing back into (|22p gives the quadratic equation 

/?2 - 3/? + 1 = , (24) 
whose only solution leading to a positive ymin is 

(3 = l + ip, V'^^^. (25) 
ip is the golden ratio. Using the recurrence property of the golden ratio 

(^"+1 =^^ + (^"-1 , (26) 

we arrive at ymin = 293, and therefore at 

\/5 — 1 

Mk = (pLmax Lmax = '^Mk , <^ = , (27) 

where $ is the golden ratio conjugate. Thus, quite strikingly, the two null circular orbits merge 
when the ratio between the Komar mass of either black hole and the semi-distance between them 
is the golden ratio! This merging happens for Weyl radial coordinate 

p = y/2 + ifMK . (28) 
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Figure 2: Left: Effective potential for i) L = 0, which has two extrema at pcNOii^^ = 0) = and 
PCNO'zi^ = 0) = ^J?iMadm] ii) L = Mx/3, for which the two extrema have moved closer to one 
another; iii) L = Lmax = = Mk/^, for which the two extrema have merged. We have fixed 

Mk = 1 in the plots. Right: Radial coordinate p of the photon circular orbits as a function of L. 



3.1.2 Circular time-like orbits 

The significance of the circular photon orbits described in the last subsection can be seen by 
considering circular time-like orbits. From (jl6p these are obtained by extremising the potential 



Vip) 



,2(7(p,0) _^ j2 



e4C/(p,0) 



(29) 



J2 

+ — 
P^ 



\ Madm + L + ^J{Madm + Ly + py \Madm + L+ ^/{Madm + Lf + p^ j 

where J is now the angular momentum per unit mass. The extrema of this potential are given by 

1 



f{p,MADM,L) 



2 + iJp)-^ y^iMADM + Ly + p2 + Madm + L ^/L^ + p^ - L 



(30) 



In Fig. O we display the solution of Eq. (j3U|) for p, in terms of J, for fixed values of L. One 
observes that for L < Lmax = l/v' there is a forbidden hand for time-like circular geodesies, which 
is the region between the two photon circular orbits. This behaviour should be contrasted with 
that exhibited in Fig. [T]for the Newtonian case and is illustrated in Fig. [H 
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Figure 3: Radial coordinate p of the circular time-like orbit as a function of the angular momentum 
per unit mass J, for various values of L and fixing Mk = 1- One observes that there is a forbidden 
region p G [pcNOi-, PCN02\^ ^ L^ax, where pcNOi are the radii of the two circular null orbits. 

The time-like orbits approach the photon orbits for J — >^ oo, from either side. 



3.2 Chcirged case 

Let us start by recalling that in a Reissner-Nordstrom black hole, the circular null orbits are found 
as extrema of the potential 

+ (31) 

in standard Schwarzschild-type coordinates, where q is the charge to mass ratio of the black hole. 
Extremising this potential gives the location of the circular null orbit outside the horizon 

rcN02 = ^ + y ( ^ 1 - ^r^ADM ' (32) 

whereas the second null orbit, with constant radial coordinate, relevant for our analysis is located 
at the horizon 

rcATOi = Madm + \]Mljjj^-q^Ml^j^ . (33) 

Thus, the forbidden band for circular time-like orbits exists for all black holes in the Reissner- 
Nordstrom family. In the extremal case, this band acquires the special property of having infinite 
area, since the horizon, at r = Madm = i^CNOn is at an infinite proper radial distance from any 
point with radial coordinate r > Madm, which is the case for r = 2Madm = 'rcN02- 
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Figure 4: Illustration of the behaviour of circular photon orbits (dotted lines) in the double- 
Schwarzschild solution. For L = 0, the two circular photon orbits sit at r = 3Madm (in 
Schwarzschild coordinates) and at the horizon (top). As L increases these two photon orbits 
approach one another (middle) and merge as L = is approached (bottom). The two photon 

orbits delimit the forbidden annular region for circular time- like orbits (coloured region). 

In order to consider two charged Reissner-Nordstrom black holes we shall now discuss charged 
Weyl solutions. The background fields are 



The electrovacuum Einstein-Maxwell equations 

1 



D^F^"" = 



-Xip,z)dt. (34) 



(35) 



give rise to the set of equations 

A^^U = e-^^ [{dpxf + [dzxf] , A^sx = 2[dpxdpU + d,xdzU] , (36) 

and 

d,k = 2pdpUd,U - 2pe-^^dpxdzX , dpk = p[[dpUf - [d,Uf] - pe-^^[{dpx? - {d,x?] ■ (37) 

This set of equations is more involved than eq. ([I])-®, in particular because there is no longer a 
linear equation to be solved. However, if one takes 



e^^'^P^'^ = l--x{p,z)+x{p,zf 



(38) 
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where g is a constant, and 



(39) 



the Einstein-Maxweh equations for (p4|) reduce to exactly the eq. (jl])-© [TO]. Thus we take U{p, z) 
and k{p,z) to be the same as in the uncharged case, ([6]) and ([9]). The background (134]) wih then 
describe a double Reissner-Nordstrom solution, wherein both black holes have the same charge to 
mass ratio q. 

The study of circular null geodesies will follow that of the uncharged case, mutatis mutandis; 
they will be given by the extrema of the potential 

4 



Vip) 



o4C/{p,0) 



=4(7(p,0) 



p^ p^ 

where e^^'''^^ may be read off from ()17p . Such extrema obey 

1 Q 



(40) 



f{p,MADM,L) 



where 



Q 



2 g + e-2W) ' 



(41) 
(42) 
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Figure 5: Radial coordinate of the photon circular orbits, R = p\/l — q^, in terms of the distance 
L for various values of the charge Q defined by (j42p . 

The result for the radii of the photon circular orbits in terms of the distance L for various values 
of the charge is shown in Fig. [5l where we used the coordinate R = p\/l — q"^. The first feature 
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we wish to emphasise is that, as L is increased for fixed charge, the two photon orbits approach 
one another, just as for the uncharged case, and coalesce for a maximum value of L. Thus, the 
forbidden band for time-like circular orbits exists for all possible values of the charge, for sufficiently 
small L. The second feature we want to mention is that, fixing L, Mk and increasing the charge, 
the forbidden band has larger area. For L = this area diverges as extremality is approached, as 
mentioned above. In Fig. [6] we exhibit an example with L 7^ 0. 

In the extremal case, for which the solution is completely regular on and outside the event 
horizon, the analysis may also be done using the Majumdar-Papapetrou |lllll2j form of the solution: 



ds^ = -tt/^ + H{p, ^)\dp' + p'dcP' + dz') , A = --^ , (43) 
H{p,zy H{p,z) 

where 

mp.z) = l+ , + , . (44) 

In this extremal case, circular null orbits exist in the z = sub-manifold as long as the equation 

(p2 ^ ^2^3/2 ^ 2Mx(p^ - L^) , (45) 

has solutions for real and positive p, which for fixed Mk is possible when L obeys 

/2\ 3/2 

L<Lmax=i-j Mk . (46) 



4 Interpretation: Optical Geometry 

A D dimensional static geometry with SO{D — 1) isometry group may always be expressed by the 
line element 

ds" = gtt{R)dt^ + gRR{R)dR^ + gee{R)dVLD-2 , (47) 

where dVtD-2 is the line element on the (Z) — 2)-sphere. The optical geometry is the effective spatial 
geometry seen by light rays: 

,2 9rr{R) ,n2 , 9ee{R) /.^^ 

defined where gu{R) < 0. 

For the spacetime (|47p . the radial equation of motion of a particle with mass, energy and angular 
momentum m,E,J, respectively is 

\gu{R)\9RR{R)R^ = E^-\gu{R)\(m^ + -^) . (49) 

Circular null orbits are therefore determined by the extrema of the potential 

geeiR) 



which is both the coefficient of the angular momentum term in ()49p and the inverse of the proper 
radius squared of the sphere line element in the optical geometry. This potential pushes the test 
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Figure 6: ^Jgpp as a function of p, for fixed Mk = 1, L = 0.5, and two values of the charge. The 
vertical lines correspond to the radii of the null circular orbits, and the shaded region corresponds 
to the area of the forbidden band for time-like circular orbits divided by 27r. In the Q = 0.99 case 
the exterior circular photon orbit is outside the p range. Observe that the area increases with the 
charge. 



particle in the direction of increasing proper size of the spheres in the optical geometry. We dub this 
direction as outwards. If it coincides with our nai've notion of "outwards" , the angular momentum 
term may still be interpreted as a centrifugal term, as in flat space. Time-like circular orbits are 
only possible in a region where this term and the mass term in (j49p originate forces in opposite 
directions. 

For a single Schwarzschild black hole of mass M, circular time-like orbits are possible for r > 3M 
in Schwarzschild coordinates. Thus they are forbidden, outside the horizon, for 2M < r < 3M. In 
this case, the optical geometry is 

rf^optical = (1 _ 2M/r)2 + l-2M/r'^^^ ' ^^^^ 

The proper size (area) of the 2-sphere is therefore A = At: /V{r), where the potential V{r) is given 
by (jl9p . The area increases with r for r > 3M, but decreases for 2M < r < 3M - Fig. [7| (left). 
Thus the "centrifugal force" points towards the black hole, for 2M < r < 3M. Since the mass term 
is always attractive towards the black hole, there can be no circular time-like orbits in the region 
2M <r < 3M. 

For the double Schwarzschild solution the optical geometry in the symmetry plane is 

<ticai = e-'^^^'''> le^'^^^'^Up' + p'dc^'] . (52) 
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Figure 7: Left (right): proper size of the 2-sphere (circle) in the optical geometry of the single 
(double) Schwarzschild solution. Defining "outwards" as the direction in which this size increases, 
then "outwards" corresponds to decreasing the radial coordinate for r < SMadm (an annulus- 
like region, when < L < $). In the single Schwarzschild case, the optical metric is singular at 
r = 2M; for r < 2M, the proper size is negative since V{r) < 0. In the double-Schwarzschild case, 
the annulus-like region lies in between the two extrema of the curves. We have set the individual 
black hole masses to unity. 



The proper size (perimeter) of the cj) circle is therefore given by £ = 2-K/y^V{p), where the potential 
V{p) is given by ()17p . The forbidden band for circular time-like orbits is the region wherein the 
radial coordinate decreases "outwards" , and therefore, the "centrifugal force" becomes directed to 
the symmetry axis - Fig. [7] (right) - preventing again the existence of time-like circular orbits. 



5 Multiple black holes 

We shall now consider a Weyl solution with N uncharged black holes. Then ([6]) is replaced by 

(53) 



^2U{p,z) _ (-^1 ~ Cl){R3 - Cs) {R2N-I - C2N-1] 
' ~ (i?2-C2)(i?4-C4) {R2N-C2N) 



This function U{p, z) describes the Newtonian potential of N rods of infinitesimal width and mass 
density g = 1/2, located at p = and with z coordinate in the intervals, respectively, 

[01,02], [03,04], ••• , [a2N-i,a2N] ■ (54) 

Imposing 

On = -a2N+i-n , n = 1 . . . 2A^ , (55) 
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the geometry still admits a discrete symmetry, of which z = is a fixed point set. Then, the 
circular null orbits in this totally geodesic sub-manifold are obtained by extremising the potential 

2 



V{p) 



„4C/(p,0) 



1 



p2 + a2jv_i + a2Ar-l 



p2 + 02JV + a2Ar 



which yields the condition 



2N 



E 



We set 



0'2n — a2n-l = 2Mr- 



(-ira« 

1 ^/p^ + al 



«2ra+l — — 2L 



n=l...N , 



(56) 



(57) 



(58) 



such that the coordinate distance between each two black holes is 2L and their Komar mass is Mk- 
In Fig. Owe have plotted the radial coordinate of circular null orbits versus L, for various numbers 
of black holes. For even, the behaviour is quite similar to that discussed for N = 2. Moreover 



lim Lr, 



Mk . 



(59) 



For N odd, the behaviour is quite similar to that of a single Schwarzschild black hole, due to the 
existence of a black hole in the symmetry plane. 




Figure 8: Radial coordinate of circular null geodesies for an array of N black holes as a function of 
L. Left (right): N odd (even). 



6 Final Remarks 

In this paper we have studied a relativistic version of Euler's 3-body problem: the motion of 
test particles in a Z2 invariant Weyl solution describing two (or more) black holes. We have 
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considered both uncharged and charged black holes, within Einstein-Maxwell theory. Generically, 
these solutions have conical singularities on the symmetry axis. In the extremal case, however, they 
become a special case of the Majumdar-Papapetrou spacetimes, and are completely regular on and 
outside the event horizon [13\ [T^ . 

For sufficiently small distance between the two black holes, there are two circular null geodesies 
in the symmetry plane corresponding to the set of fixed points of the Z2 symmetry. The interior 
one is stable (against perturbations in the symmetry plane), while the exterior one is unstable. 
These orbits delimit a forbidden region for circular time- like geodesies. In the optical geometry, 
the proper size of the 50(2) isometry group orbits grows towards the symmetry axis, rather than 
away from it, in this region. This unusual property is correlated with a physical effect that may 
be, heuristically, described as the centrifugal force becoming directed towards the symmetry axis. 
Thus, it cannot balance the gravitational attraction, preventing the existence of time-like circular 
orbits. 

The existence of a stable circular null geodesic is associated with a behaviour which may be 
described as an 'optical fibre' or 'wave guide' like geometry. Photons with sufficiently low energy 
will be caught in a potential well in the radial direction, as may be observed in Fig. [2j Thus, as 
they move in the (f) direction, photons will oscillate between a maximum and a minimum radii. The 
wave guide like geometry may allow light waves to travel around the (p direction without decreasing 
their amplitude significantlyo 

As the distance between the two black holes is increased, the two null circular orbits approach 
one another, eventually coalescing. For the uncharged case this happens when the ratio between 
the Komar mass of each black hole and the semi-distance between the black holes is the golden 
ratio. 

One may ask whether the geodesic equations in the full spacetime, rather than just the symmetry 
plane, are Liouville integrable. In [T7j a similar question was studied for extremal black holes in 
Einstein-Maxwell-dilaton theory, for arbitrary dilaton coupling. These authors showed that the 
geodesic motion is generically chaotic. There is only a very special case of integrable motion: 
when the dilaton coupling is that of Kaluza-Klein theory. For the generic case studied herein, of 
charged Weyl solutions that may be continuously connected to the extremal black holes of Einstein- 
Maxwell theory (with zero dilaton coupling therefore), it seems quite unlikely that the motion 
may be integrable. Observe however, that the uncharged double Schwarzschild solution may also 
be continuously connected to the extremal black holes of Kaluza-Klein theory. Thus, it would 
be interesting to understand where, along this trajectory in the space of solutions, integrability 
appears. 
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